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Data from the experiment on the reaction at 23 AMeV are analyzed with an aim 

to find signatures of exotic nuclear configurations such as toroid-shaped objects. The experimental 
data are compared with predictions of the ETNA code dedicated to look for such configurations and 
with the QMD model. A novel criterion of selecting events possibly resulting from the formation of 
exotic freeze-out configurations, ’’the efficiency factor”, is tested. Comparison between experimental 
data and model predictions may indicate for the formation of fiat/toroidal nuclear systems. 

PACS numbers: 25.70.Gh, 25.70.Pq 


I. INTRODUCTION 

The search for exotic nuclear configurations was in¬ 
spired by J.A.Wheeler [ij. His idea was investigated by 
many authors who studied the stability of exotic nuclear 
shapes (see e.g. [l-Q)- Theoretical investigations have 
shown that very exotic extra superheavy nuclei can be 
reached only if non-compact shapes are taken into ac¬ 
count. Calculations for bubble structures showed that 
such nuclei can be stable for Z > 240 and N > 500 
(see e.g. [B-I3|)- Recently it was found that for nuclei 
with Z > 140 th e g lobal energy minimum corresponds to 
toroidal shapes Q- In contrast to bubble nuclei, the 
synthesis of toroidal nuclei is experimentally available in 
collisions between stable isotopes. 

To address this issue simulations were performed for 
Au + Au collisions in a wide range of incident energies 
using the BUU code ©[HI- These calculations indicate 
that the threshold energy for the formation of toroidal 
nuclear shapes is located around 23 AMeV. 

Also Improved Quantum Molecular Dynamics Model 
calculations performed for U J- U collisions have shown 
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a possible formation of toroidal freeze-out configurations 
above a specific collision energy for this heavy system 
[l^ . Such toroidal-shape complex can be also created 
in macroscale in binary droplet collisions above some 
threshold velocity [l^ . 

A number of observables were suggested as the sig¬ 
natures of noncompact freeze-out configurations. These 
were: 

• Larger number of intermediate mass fragments 
should be observed than would be expected for the 
decay of a compact object; 

• Enhanced similarity in the charge and size of the 
fragments should also be observed; 

• Suppressed sphericity in the emission of fragments 
should be visible. 

The simulations of decay process of different break up 
configurations using the ETNA code were performed to 
study the ability of the CHIMERA detector [3, [I^ for 
recognition of non-compact configurations. Analysis of 
different observables have shown that a quantity named 
“the efficiency factor” of events with 5 heavy fragments 
can be used as a criterion of selecting events possibly 
resulting from formation of toroidal configurations |I6l . 

El- 
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The CHIMERA collaboration has carried out an exper¬ 
iment on the reaction at 23 AMeV beam 

energy, focused on two independent goals, first on the ex¬ 
tension of the earlier study at lower energy of 15 AMeV, 
in which a new reaction mechanism of violent breakup of 
non-fusing -I- system into 3 and/or 4 mas¬ 

sive fragments was observed jl^, [l^, [l^, and second, 
on the search of exotic nuclear configurations such as 
toroidal shapes. Some preliminary results of the former 
project have been published in f^ . 

In this work we report results of our analysis focused 
on the question of exotic configurations involving the 
breakup of the -I- Au system into 5 or more 

fragments. The experimental data are compared with 
model predictions. Conclusions regarding the shape of 
the freeze-out configuration are drawn. 

This paper is organized as follows. In Sec. 2 we present 
the experiment and data calibration procedures. General 
characteristics of experimental data are shown in Sec. 3. 
The dedicated observables are discussed in Sec. 4. The 
conclusions are presented in Sec. 5. 


II. EXPERIMENT AND DATA CALIBRATION 
PROCEDURE 

The experiment for the -|- reaction at 23 

AMeV was performed at INFN-LNS Superconducting 
Cyclotron of Catania. During the experiment two gold 
targets were used: 164 and 396 ^g/cm?. The thinner tar¬ 
get was used in calibration measurements and the thicker 
one in the production runs. Reaction products were de¬ 
tected with the CHIMERA multidetector [3, [3| that is 
constituted by 1192 telescopes arranged in 35 rings in full 
27r azimuthal symmetry around the beam axis, covering 
the polar angle between 1° and 176°. A single detec¬ 
tion cell is constituted of a planar n-type silicon detector 
(= 300 /i thickness) followed by a Cs(Tl) scintilator of 
thickness varing from 12 cm at forward angles to 3 cm at 
backward angles. 

The collected data were calibrated using a set of ded¬ 
icated programs developed at INFN-LNS. Energy cali¬ 
bration of Si detectors was performed using ion beams, 
delivered both by the tandem and the cyclotron. Data 
for the following systems were used: (i) the elastic scat¬ 
tering data for -|- Au at 60 and 80 MeV, -|- Au 
at 142 MeV, Au -k An at 170 MeV and 23 AMeV; (ii) 
recoil peak for Au -|-^^ C at 170 MeV; and (ii) fission 
fragments from Au C reaction at 23 AMeV. Unlike 
in the analysis of ternary breakup reactions [3| in which 
the pulse-height defect was calculated with the formula 
of Ref. [13] , in the present analysis the pulse-height de¬ 
fect in silicon detectors was calculated using the same 
procedure as described in Ref. [Isj ]. 

In order to identify fragments two methods were ap¬ 
plied: (i) the AE — E technique for fragments punch¬ 
ing through the silicon detectors; (ii) the time of flight 
(TOE) method for the class of fragments stopped in Si 



FIG. 1: AE — E spectrum(upper panel) and the 
corresponding Z spectrum (bottom panel) for fragments 
detected in telescope placed at 0 = 5.2° for Au -|- Au 
reaction at 23 AMeV. 

detectors. In Fig. [T] (upper part) an example of AE — E 
plot is shown for a detector belonging to the 3-th inter¬ 
nal ring, at a polar angle 9 = 5.2°. The Z distribution 
of fragments identified by the AE — E method for the 
same detector is presented also in Fig. [T] (bottom panel). 
We can see that in the Z spectrum a good charge identi¬ 
fication can be observed up to Z = 42. At this angular 
range one observes a broad maximum of the charge dis¬ 
tribution located at charges Z = 30 — 40 corresponding 
to Au fission fragments. One observes also a substantial 
contribution of lighter fragments. In order to estimate a 
missing information on mass of the fragments identified 
in charge by AE — E method, the EPAX formula mil 
was used. 

The mass of fragments stopped in Si detector is de¬ 
termined by TOF method. The start signal was given 
by 30% Constant Fraction Discriminator acting on time 
signal generated by the silicon detector, while the stop 
signal was given by delayed Reference Signal delivered by 
cyclotron. Examples of AE — TOE spectra are presented 
in Fig. [5J In this case mass values are calculated using 
the formula: 

m = 2E-(to- (1) 

where R is the distance between the target and a given 
detector and the to is a time offset of the measured time 

t. 

A crucial problem in the calibration of TOF measure¬ 
ments for the CHIMERA multidetector is evaluation of 
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FIG. 2: (color online) The Ai? — TOF distributions for 
detectors 404 (upper panel) and 653 (bottom panel) 
located at 0 = 17° and 28.5°, respectively. Lines 
presented of both figures correspond to positions of 
masses as indicated. 


to offset that must be determined for each detector in¬ 
dividually. Moreover, to depends on mass, charge and 
kinetic energy of the detected fragment. The to values 
for well identified light fragments and Au-like nuclei frag¬ 
ments located at the left-hand-side edge of the AE—TOF 
distribution (see Fig. [2|) are presented by color symbols 
in Fig. [21 For relatively light fragments a well tested 
parametrization of to for CHIMERA detectors was pro¬ 
posed [^ . To calibrate tg for medium and heavy frag¬ 
ments in -b experiments a new calibration 

method based on a functional dependence of to on mass, 
charge and pulse-height-defect dependent kinetic energy 
was developed (see e.g. [2^) and applied in analysis of 
the ternary breakup experiment 211. 

In the present analysis we use another method of the 
parametrization of to offset: 


I to^sat to^sat ^ At 

I to,sat At to,sat ^ At 




Energy (MeV) 

FIG. 3: (color online) The to dependence on incident 
energy and particle mass for telescopes 404 and 653 
located at 6* = 17° and 28.5°, respectively. Color 
symbols represent the to values for identified light 
fragments and Au-like fragments. Solid lines represent 
the htting results using the formula. The dashed line 
indicate the to,sat value. 


The values of the parameters A, B, and e 

were determined by fitting Eq. (2) to selected points for 
well identified light fragments and Au-like nuclei frag¬ 
ments (see Fig. I2D. In this procedure fragment energies 
were corrected for pulse-height defect and their Z values 
were estimated using the EPAX formula [22-24]. Fig. [3] 
presents the to calibration lines calculated for selected 
mass values. The calculated positions in the AE — TOF 
distributions for some selected mass numbers are shown 
for two detectors in Fig. [3] as solid color lines. 

In Fig. m two dimensional mass versus kinetic energy 
distributions are shown for telescopes located at two an¬ 
gular regions. For 3° < 0 < 10° region (upper panel) 
the distribution extends from small masses seen at low 
energies up to the Au elastic peak. For 20° <9 < 29° re¬ 
gion (bottom panel) the particles with masses up to 200 
a.m.u. are observed at relatively low kinetic energies. 


At = B-A(l- expij ■ m)) ■ (^)(«-<5-”^) 

EpT 

hpT 


(3) 


where to,sat is determined for particles punching 
through the silicon detector. The Epr is the highest 
energy deposited by particles with mass m. 


III. THE GENERAL CHARACTERISTICS OF 
EXPERIMENTAL DATA 


In Fig.|5]two dimensional distribution mass versus par¬ 
allel velocity of identified fragments is shown. Location 
of quasielastic Au peak is visible at mass around 200 and 
velocities close to the beam velocity {vp = 6.67cm/ns). 
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FIG. 4: (color online) The correlation between mass 
and energy for fragments observed in telescopes located 
at 3° < 0 < 10° (upper panel), and 20° < 6» < 29° 
(bottom panel). 


Peak corresponding to An recoil fragments can be found 
at velocities close to zero. Here one can observe an un¬ 
derestimation of mass value for these fragments due to 
the imperfection of used to parametrizations (see Eq. 2). 
At velocities between these two limits fragments origi¬ 
nating from fission of the Au-like nuclei are located. One 
can also identify a separated region located at low masses 
and velocity close to center of mass velocity. This region 
correspond to the intermediate velocity source. 

For the registered events we have constructed the plot 
presenting the dependence between the total charge of 
identified fragments, Ztot^ versus total parallel momen¬ 
tum of those fragments normalized to the beam momen¬ 
tum, Pw^tot/Pproj (see Fig.H]). 

One can distinguish different regions on this plot. In 
the region of low values of total collected charge and low 
paralled momentum one observes the ridge corresponding 
to the badly detected events. In the region of total paral¬ 
lel momentum close to 1 and total collected charge close 
to the charge of projectile one observes the maximum 
corresponding to deep inelastic collisions where the tar¬ 
get like fragment remains undetected. Region where the 
total detected charge is close to total charge of the system 
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FIG. 5: (color online) The correlation between mass of 
identified fragments versus parallel velocity of those 
fragments. 
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FIG. 6: (color online) The correlation between total 
charge of identified fragments versus total parallel 
momentum of those fragments normalized to the beam 
momentum. 


and the total parallel linear momentun is close to linear 
momentum of the projectile can be called as region of well 
reconstructed events. In our present analysis this region 
is selected imposing the conditions: 120 < Ztot < 180 
and 0.8 < p\\^tot/Pproj < 1-1. The number of events full¬ 
filing these conditions is equal 5.9 * 10®. 

For this class of well reconstructed events in the Au 
+ Au reaction the multiplicity distribution of fragments 
with charge Zfrag > 3 and Zfrag > 10 are presented 
in Fig. [3 One can notice here that the number of 
events with five or more fragments corresponding to 
above charge thresholds is equal about 116 000 and 6 000, 
respectively. 


IV. DATA ANALYSIS 

Based on results of Ref. [l^ , [l^ and , in analysis 
of ternary and quaternary events one expects the ob- 
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FIG. 7: (color online) Multiplicity distributions of 
fragments with Zfrag > 3 (red histogram) and 
Zfrag > 10 (blue histogram), respectively. 


servation of binary deep-inelastic collisions followed by 
breakup of one or both primary reaction products. Re¬ 
sults of analysis of this particular class of reactions can 
be found in j^ . 

For the class of events with five fragments one can con¬ 
sider at least two mechanisms responsible for the presence 
of the fifth heavy fragment: (i) creation of the fragment 
in the interaction region (intermediate velocity source) 
for peripheral collisions or (ii) the multifragmentation of 
the composite nuclear system formed in central collisions. 

In order to investigate the reaction scenario responsible 
for events with five and more fragmets we have compared 
experimental data with ETNA and QMD model predi- 
cions. The ETNA model can simulate the decay of nu¬ 
clear system assuming compact and noncompact freeze 
out configurations |l6l |. In this model three freeze out 
configurations are considered: (i) ball geometry with vol¬ 
ume 3 and 8 times greater than normal nuclear volume 
Vo (fragments uniformly distributed inside the sphere); 
(ii) fragments distributed on the surface of the sphere 
mentioned above (bubble configuration); (iii) fragments 
distributed on the ring with diameter 12 fm and 15 fm 
(toroidal configuration). In this model we consider events 
corresponding to central collisions only (0-3 fm impact 
parameter range). 

In order to simulate the contribution from noncen¬ 
tral collisions the QMD model [2^ calculations were per¬ 
formed in the full impact parameter range 0-12 fm. In 
our analysis the QMD code developed by Lukasik et al. 

was used. This code takes into account: (i) protons 
and neutrons (in standard QMD each nucleon has an 
effective Z/A charge); (ii) momentum dependent Pauli 
potential (Skyrme + Coulomb -I- Symmetry -|- Surface 
-I- Pauli) is used instead of the Yukawa potential; (iii) 
initial nuclei in their real ground states with minimum 
energy (thanks to the Pauli potential) are prepared,(iv) 
strict angular momentum conservation in collisions is ap¬ 
plied; (v) to simulate the Pauli blocking: in Lukasik code 
the overlap of 6 -dimensional Gaussians is used, unlike in 
standard QMD where overlap of appropriate spheres in 



FIG. 8 : (color online) The coplanarity vs sphericity 
distributions for Ball 8 Vo, Toroid 15 fm, QMD and 
experimental data. 


configuration and momentum space is included. 


A. The shape sensitive observables 

In our analysis several observables sensitive to the 
freeze-out break-up configuration are investigated. As 
a first, we consider the shape of events in the momentum 
space M- The diagonalization of the momentum tensor 
gives three eigenvalues Ai and three eigenvectors ef. The 
sphericity and coplanarity variables are defined as: 


s = 1.5(1 -Ai),c=^(A 2 -As), (4) 

where Ai > A 2 > A 3 are normalized to their sum. 

In the coplanarity vs sphericity plane all events are lo¬ 
cated inside a triangle defined by points ( 0 , 0 ), (f,-^), 
and (1,0). In Fig. [ 8 ] the (s,c) distribution for experimen¬ 
tal data is compared to the ETNA model predictions 
for Ball 8 Vb, Toroid 15 fm freezeout decay configurations 
and with QMD predictions. In the case of ball geom¬ 
etry the maximum of the corresponding distribution is 
located in the centre of the triangle. For toroidal con¬ 
figuration the distribution is located closer to the line 
(0,0), (|, ^). One can see that the experimental distri¬ 
bution looks very similar to QMD distribution which is 
dominated by noncentral collisions contribution. 

In order to reduce noncentral contribution we have in¬ 
vestigated for the QMD model predictions the depen¬ 
dence between flow angle, 0fiow, and impact parameter 
(see Fig. (panel a), where 9flow is the angle between 
beam axis and the eigenvector et for the largest eigen¬ 
value Al. One can see on this plot that most noncentral 
events are located at small 9flow angles. The similar 
dependence is observed for experimental data between 
9flow and total transverse momentum, ptrans, used as 
impact parameter estimator (see Fig. [S] (panel c). We 
decide to reduce contribution of noncentral events both 
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FIG. 9: (color online) The impact parameter vs 9fiow 
dependence for QMD model predictions (panel a) and 
total transeverse momentum, ptrans, vs 9flow 
dependence for experimental data (panel c). The 
impact parameter vs 9piane dependence for QMD model 
predictions (panel b) and total transeverse momentum, 
Ptrans, VS 9plane dependence for experimental data 
(panel d). 


for experimental data and model predictions by using the 
condition 9flow > 20°. 

The (5, and observables as most sensitive to the 
shape of freeze out configurations were selected [T^. The 
5 variable is related to sphericity and coplanarity vari¬ 
ables . The 5 variable measures the distance between 
a given point of the (s,c) distribution and the line de¬ 
fined by points (0,0), (|, -^). In the Fig. [TO] (left panels) 
the 5 distributions are presented for experimental data, 
ETNA model predictions for considered freeze-out ge¬ 
ometries and QMD predictions. One can see here that 
the 5 distribution for experimental data is similar to that 
corresponding QMD predictions. The biggest difference 
can be observed with the distribution for Ball 8Vb con¬ 
figuration. 

The variable used in our analysis gives a measure 
of the event flatness in the velocity space. For each event 
we are establishing the plane in the velocity space. The 
parameters of this plane are selected in the way that 
the sum of squares of distances between the plane and 
the endpoints {vx,i, Vyf, Vz,i) of velocity vectors reach the 
minimum value. This last quantity is called the A^ pa¬ 
rameter and is defined as: 

Nfr 

A^ = min[J2(.dUA,B,C,D)], (5) 

i=l 

where: 

^ _ \A * Vx^i B * Vyf C * Vzf D\ 

and parameters A, B, C, and D are the plane param¬ 
eters. The plane parameters and the velocities of frag¬ 
ments are in the velocity of light units. 



FIG. 10: (color online) In the upper left panel the S 
distributions are presented for experimental data. Ball 
8Vb) Bubble 8Vb freeze-out geometries and QMD 
predictions. In the bottom left panel the experimental 
distribution is compared with predictions for Toroid 12 
fm and Toroid 15 fm configurations. In the right panels 
the A^ distributions for experimental data and model 
predictions are shown. All the distributions presented 
here are constructed using the condition Zfrag > 10 and 
9fiow > 20 °. 


The A^ distributions are shown in Fig. [TU] (right pan¬ 
els) for data and model predictions. One can see here that 
for A^ variable the biggest difference between experimen¬ 
tal distribution and model predictions is observed for the 
Ball 8Vb, and Bubble 8Vb conhgurations. In contrast to 
that, the experimental data seem to be more consistent 
with the simulations assuming toroidal freeze-out config¬ 
urations. 

In relation with A^ parameter one can define an angle, 
9 plane, between the beam direction and vector normal 
to the plane defined by parameters A, B, C, and D. For 
events corresponding to noncentral collisions, where most 
of reaction products are located in the reaction plane, 
9piane should be close to 90°. This behavior is illustrated 
in Fig. ini (panel b) for QMD model predictions, where 
most of noncentral events are located in the reaction 
plane. The similar dependence is observed for experimen¬ 
tal data between 9 flow and total transverse momentum, 
Ptrans, used as impact parameter estimator (see panel d). 

The dependence between 9 plane and 9 flow for Ball 8Vb, 
Toroid 15 fm, QMD and experimental data is presented 
in Fig. [TT] One observe here that for experimental data 
most of events is located in the region selected by condi¬ 
tions 9flow < 20° and 9piane > 75°. The same behavior is 
observed in the case of QMD calculations. These observa¬ 
tions indicate that such events correspond to noncentral 
collisions. For the Ball 8Vb conhguration one observes the 
correlation between 9flow and 9piane angles. For toroidal 
conhguration the correlation between these angles is even 
stronger. Most of these events is located in the region de- 
hned by conditions 9flow > 20° and 9piane < 75°. 

Following the method proposed in Ref. [l^ we select 
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FIG. 11: (color online) The dependence between 9plane 
and 9flow for Ball 8Vb (panel a), Toroid 15 fm (panel 
b), QMD (panel c) and experimental data (panel d). 


events corresponding to a toroidal shape by the set of 
conditions: 

< 0.001 and S < 0.05. (7) 

As an efficiency measure of the above conditions we 
take ratio of number of events fulfilling the selection con¬ 
ditions to the number of events with five and more heavy 
fragments. Hereafter, this ratio is called the efficiency 
factor (EF). 

The results of this procedure are presented in the Fig. 
M for different regions of 9flow and 9piane angles. As 
one can see the EF is very low for spherical freeze-out 
configurations with respect to the corresponding values 
for toroidal configurations. 

For QMD calculations the value of the efficiency factor 
is strongly dependent on the 9piane range. The condition 
9piane < 75° reduccs the number of flat noncentral events 
mostly located in the reaction plane. For events selected 
additionally by the condition 9flow < 20° the EF drops 
to zero. 

For experimental data the value of the efficiency factor 
is about 50% for events located in the reaction plane 
{9piane > 75°) and is reduced by factor of 2 for events 
perpendicular to the reaction plane. These values are 
weakly dependent on the 9 flow angle range. 

One observes that the values of the EF for experimen¬ 
tal data are much larger than the correspondig predic¬ 
tions for QMD model. The biggest difference is observed 
for events located outside the reaction plane (9piane < 
75°) at small 9 flow angles. 

In order to investigate a possible formation of toroidal 
configurations in our analysis we selected the region 
where according to ETNA predictions the toroidal config¬ 
uration is most pronounced in the 9flow and 9plane plane 
{9piane < 75° and 9 flow > 20°). In Table Uthe efficiency 
factor values are given for experimental data and model 
predictions. The efficiency factor values are shown for 
four threshold values of the fragment charge. 


H Gplane < 75° 



FIG. 12: (color online) The EF values for different 
windows of 9plane and 9flow- The presented results were 
sorted using the condition Zfmg >10. 


Efficiency factor (%) 


Configuration 

Zfrag ^ 3 

Zfrag > 10 

Zfrag > 15 

Zfrag > 20 

Ball 3Vo 

3.3 ± 0.2 

3.5 ± 0.2 

3.5 ± 0.2 

3.5 ± 0.2 

Bubble 3Vb 

2.4 ± 0.2 

2.6 ± 0.2 

2.7 ± 0.2 

2.7 ± 0.2 

Ball 8Vo 

3.2 ± 0.2 

3.5 ± 0.2 

3.5 ± 0.2 

3.5 ± 0.2 

Bubble 8Vb 

3.9 ± 0.2 

4.6 ± 0.2 

4.7 ± 0.2 

4.7 ± 0.2 

Toroid 12 fm 

29.7 ± 0.6 

31.6 ± 0.6 

31.8 ± 0.6 

31.9 ± 0.6 

Toroid 15 fm 

25.2 ± 0.5 

27.5 ± 0.5 

27.7 ± 0.5 

27.8 ± 0.5 

QMD 

13.7 ± 3.4 

8.2 ± 4.7 

6.3 ± 5.5 

N/A 

data 

27.1 ± 0.7 

26.2 ± 2.5 

26.2 ± 4.8 

21.1 ± 8.0 


TABLE I: The efficiency factor at incident energy 23 AMeV for 
four threshold values of the fragment charge for events selected by 
conditions 6/iou, > 20° and Opiane < 75°. 


From Table U we notice that the EF values for ex¬ 
perimental data are very close to the model predictions 
for toroidal configurations. This observation may be one 
of arguments in favor of the formation of toroidal/flat 
freeze-out configuration created in the Au + Au colli¬ 
sions at 23 AMeV. 



































FIG. 13: (color online) The distributions of standard 
deviation of the fragment mass for non-flat events (red 
lines) and flat events (green lines) for experimental data 
(solid lines) and QMD model predictions (dashed lines). 

All the distributions presented here are constructed 
using the condition Zfrag >10. 



FIG. 14: (color online) The distribution of relative 
velocities Vij of fragments pairs for non-flat events (red 
lines) and flat events (green lines) for experimental data 
(solid lines) compared with: Ball 8Vb and Toroid 15 fm 
(upper panel) and QMD model predictions (bottom 
panel). All the distributions presented here are 
constructed using the condition Zfrag > 10. 


B. Other observables 

In order to get additional evidence to support the hy¬ 
pothesis that toroidal objects are created the behaviour 
of other observables was investigated. We consider here 
for each event separately: (i) standard deviation of frag¬ 
ment mass {o' a), (ii) relative velocities of fragments pairs 
(vij), (iii) mean velocities of fragments as a function of 
their mass. 

First we construct these observables for events selected 
by conditions 9flow > 20° and 9piane < 75°, where obser¬ 
vation of toroidal freeze-out configurations is expected. 
The distributions of these observables are generated for 
flat events selected by condition (7) (thick green his¬ 
tograms) and non-flat events (thin red histograms) se¬ 
lected by condition: 


> 0.001 and S > 0.05. (8) 

Gomparison of the oa distributions (Fig. [T51) for flat 
and non-flat events indicates that in the case of flat events 
this distribution is slightly shifted to larger values. This 
observation is in contrast with the expectation that for 
the flat events the enhanced similarity in the size of frag¬ 
ments should be visible. The corresponding distributions 
for QMD calculations are similar (dashed lines). Their 
centroids are shifted to smaller values with respect to 
experimental data. 

In Fig.fHlone observes that the distribution of relative 
velocities for flat events is shifted to smaller velocities in 
respect to non-flat events. The corresponding distribu¬ 
tions for Toroid 15 fm and Ball 8Vb ETNA model predic¬ 
tions show a similar dependence. This observation may 
indicate that the behaviour of these Vij distributions is 
insensitive to the shape of the freeze-out configuration. 

In Fig. [12] the distributions of mean velocities of frag¬ 
ments as a function of their mass for a flat and non-flat 



FIG. 15: (color online) The distributions of mean 
velocities of fragments as a function of their mass for 
non-flat events (red lines) and flat events (green lines) 
for experimental data (points with error bars) and 
QMD model predictions (dashed lines). All the 
distributions presented here are constructed using the 
condition Zfrag > 10 


events are presented. On can observe that for flat events 
velocities of fragments decrease weaker with mass com¬ 
paring to the same dependence for non-flat events. Gom¬ 
parison with same dependences presented for Pb -I- Ag 
and Pb -I- Au systems at 29 AMeV [s^ indicates that 
toroidal configurations may by created for some subclass 
of flat events. 
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Observable 

threshold 

^flow ^ 20 
Opiane ^ *^5 

& flow ^ 20 
Opiane ^ "^5 

^fiow < 20° 
Opiane 75 

^flow ^ 20 
Opiane ^ 75 

(TAiCL.m.U.) 

^frag ^ 3 
Zfrag > 10 
Zfrag ^ 15 
Zfrag > 20 
Zfrag > 25 

72.09 ±0.47 

47.01 ±1.88 

38.31 ±2.98 

31.01±6.68 

17.51±5.07 

71.09 ±0.35 

47.15 ±1.33 

38.53 ±1.24 

31.10 ±2.60 

18.95 ±5.86 

76.43±0.52 

46.41 ±2.68 

35.24 ±5.11 

25.15 ±5.15 

20.94 ±4.82 

73.33±0.13 

45.06 ±0.58 

35.58 ±0.98 

27.17 ±1.59 

18.50 ±2.23 

Vij {cm/ns) 

Zfrag > 3 
Zfrag > 10 
Zfrag ^ 15 
Zfrag > 20 
Zfrag > 25 

3.01 ±0.01 

3.13 ±0.05 

3.16 ±0.08 

3.14 ±0.24 

2.98 ±0.31 

3.17 ±0.01 

3.30±0.03 

3.27 ±0.05 

3.25±0.11 

3.28 ±0.33 

3.27 ±0.02 

3.30 ±0.08 

3.27 ±0.15 

3.24 ±0.52 

3.26 ±0.81 

3.36 ±0.01 

3.51 ±0.02 

3.49 ±0.04 

3.50 ±0.06 

3.46 ±0.13 


TABLE II: The mean values of mass standard deviation of 
the fragments, and of relative velocities Vij of fragments pairs 
for fiat events located in different windows of Ofio^ and Opiane 
angles. 


Properties of flat events in the region where obser- 
vation of toroidal freeze-out configurations is expected 
{9flow > 20° and Opiane < 75°) can be also compared 
with properties of flat events corresponding to other re¬ 
gions of 9flow and 9plane angles. Here the considered 
regions are the same as presented in Fig. [T^l The distri¬ 
butions for a A of fragments, and Vij of fragments pairs 
are presented in Fig. [12] using the condition Zfrag > 10. 
The mean values of these distributions are listed in Table 
IE We can notice here that the corresponding mean val¬ 
ues of the distribution of a a are similar for all 9 flow and 
9piane windows for a given threshold value of the frag¬ 
ment charge Zfrag- Such observation shows us that in¬ 
formation carried by a a can not be used as an indication 
of toroidal objects formation. For Vij distributions one 
observes that the mean values for class of events located 
outside the reaction plane are smaller in comparison to 
the case of events located in the reaction plane. The 
smallest mean values are seen for the region where obser¬ 
vation of toroidal freeze-out configurations are expected. 
This observation may be used as an indication that for 
events located outside the reaction plane freeze-out con¬ 
figuration is more extended in comparison with that for 
events located inside reaction plane. 


Results obtained for the considered observables sug¬ 
gest that the formation of toroidal configurations can be 
related to a fraction of flat events tilted with respect to 
the reaction plane {9piane < 75°). The probability for 
these events is much greater than the prediction of the 
QMD model. The nature of these events should be in¬ 
vestigated. 


V. SUMMARY 


We presented an analysis of events produced in Au 
-I- Au collisions at 23 AMeV. Basic information about 
data c alibration procedure were summarized. The bu lk 


0.2 
0 . 151 - 

0.1 i 

O.OS =- 


1 = 


> 75 ’ & 0 ,|„. 
< 75 ’ & 0 ,.„. 



FIG. 16: (color online) The distributions of standard 
deviation of the fragment mass (upper panel), and the 
distributions of relative velocities vtj of fragments pairs 
(bottom panel) for flat events. The red dashed lines 
corresponds to events located inside the reaction plane 
and the green solid lines correspond to events located 
outside the reaction plane. All the distributions 
presented here are constructed using the condition 
Zfrag > 10 


properties of the experimental data were shown. The ex¬ 
perimental data were compared with ETNA and QMD 
model predictions. Proximity of efficiency factor values 
for experimental data and toroidal freeze-out configura¬ 
tions may be used as an indication of the formation of an 
exotic freeze-out configuration. The juxtaposition of the 
standard deviation of fragment mass values for events lo¬ 
cated outside and inside the reaction plane are not sug¬ 
gestive of a toroidal freeze-out configuration. Compar¬ 
ison of distributions of relative velocities for event with 
different orientation in respect to reaction plane gives evi¬ 
dence that the freeze-out configuration is more extended 
for events located outside reaction plane. The behav¬ 
ior of mean velocities of fragments as a function of their 
mass for flat and non-flat events gives an indication that 
toroidal configuration may be created for some subclass 
of flat events. 

The probability of apperence of these flat events is 
much greater than the prediction of the QMD model. 
The nature of flat events tilted with respect to the reac¬ 
tion plane should be investigated. 


Assuming that the total number of detected events cor¬ 
responds to 80% of total reaction cross section, the cross 
section related to creation of flat tilted events located in 
the region where observation of toroidal freeze-out con¬ 
figurations is expected can be estimated to be equal 17fib. 
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